Unusual doping and temperature dependence of photoemission spectra from 

manganites 
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A recent, major, puzzle in the core-level photoemission spectra of doped manganites is the ob- 
servation of a 1-2 eV wide shoulder with intensity varying with temperature T as the square of the 
magnetization over a T scale of order 200K, an order of magnitude less than electronic energies. 
This is addressed and resolved here, by extending a recently proposed two electron fluid i — b model 
for these systems to include core-hole effects. The shoulder arises from a rapid redistribution of e g 
electron density, as a function of T, between the highly localized(^) and band-like (6) states. Fur- 
thermore, our theory leads to a correspondence between spectral changes due to increasing doping 
and decreasing T, as experimentally observed. 



Doped rare earth manganites, Ri^ x A x MnO^ (R = 
rare earth ion; A = alkaline earth ion), show exotic 
properties such as colossal magneto-resistance (CMR) 
close to a temperature (T) driven Ferro-metal (FM) to 
para-insulator (PI) or para-metal (PM) transition, an 
extraordinarily rich phase diagram with varying x, T 
and magnetic field (H), etc [2j. Particularly intriguing 
among their properties is the redistribution of spectral 
weight over an energy window of several eV in valence- 
band [1, [3] and Mn 2p core-level 0, @ photoemission 
spectra (PES), as a function of x and of T. The changes 
in the spectral features with doping, incompatible with 
even an approximately rigid-band description, suggests 
qualitative and drastic modifications of the underlying 
electronic structure with small changes in x. Similarly, 
changes in T affect the spectra over energy ranges 100 
times the thermal energy scale and cannot be under- 
stood in usual terms. Curiously, changes in spectra with 
increasing x or with decreasing T show striking simi- 
larities [J], suggesting a common mechanism. The ef- 
fects have been attributed 0, 0| to the unusually strong 
electron-spin and electron-lattice couplings in these com- 
pounds, but no specific theory has been proposed till 
now. 

Conventional core-level photoemission calculations, in 
terms of cluster models [7] or impurity models Q, can 
not account for a drastic renormalization of the underly- 
ing electronic structure with changing x or T. The only 
known mechanisms Q that yield T dependent changes 
in PES are Fermi-edge decoherence (FED) effects due to 
thermal excitations and Debye- Waller type effects due 
to the scattering of electrons by phonons. They lead to 
changes only on energy scales of order T and the Debye 
temperature respectively (~ 10-100 meV), and not over 



several eV as observed. Cluster calculations [f| with ad- 
justable fitting parameters, using an MnOe octahedron 
coupled to a single level at the Fermi energy, reproduce 
the multiplet features of the observed spectra; but no x or 
T dependence, of the magnitude seen in the experiments, 
can arise unless the fitting parameters are artificially al- 
lowed to vary with x and T. 

In this paper we adopt a radically different approach 
complementary to the above schemes, in that we focus 
on the x and T dependence, while neglecting details of 
multiplet structure. We employ the recently proposed 
two- fluid '£-b' model (lfj| in a dynamical mean-field the- 
ory (DMFT) [ll[ framework which successfully explains 
several hitherto poorly understood low energy properties 
of doped manganites. The active degrees of freedom in 
the manganites are the twofold degenerate e g levels, the 
t2 g core-spins of Mn, and the Jahn- Teller (JT) optical 
phonon modes of the MnOe octahedra. There are three 
strong on-site interactions, viz. the JT electron-phonon 
coupling which splits the two e g levels by an energy 2Ejt 
(~ 0.5 - 1 eV), the ferromagnetic Hund's coupling Jh be- 
tween the t2 g and e g spins (~ 2 eV) and the e g electron 
Coulomb repulsion Udd (~ 5 - 7.5 eV)[I3]; all larger than 
the e g inter-site hopping (t ~ 0.2-0.4 eV) [H The '£-V 
model (lpj is an effective low energy Hamiltonian which 
implicitly captures the crucial effects of these interactions 
and the quantum dynamics of the JT phonons. It invokes 
two types of e g electrons, one polaronic and localized (£), 
and the other band-like and mobile (b), and is given by, 
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The polaronically trapped l t species has site energy 
—Ejt, and an exponentially reduced hopping (~ 1 meV) 
which has been neglected, while the non-polaronic '6' 
species (site energy 0) has undiminished hopping t ~ 
0.2-0.4 eV . Jf is a novel ferromagnetic virtual double- 
exchange (VDE) coupling (~ 2 meV) between the core 
spins, which arises naturally in this model 0. The 
chemical potential, u, imposes the doping determined fill- 
ing constraint: J2A( n ^) + ( n b<r)) = (1 - x). 

The model is similar to the Falicov-Kimball model 
[3| , and is exactly soluble in the DMFT [llj framework, 
with the <2g core spins (Si) being approximated as classi- 
cal vectors (SCli) and Jf treated in the Curie- Weiss mean 
field approximation (loj . The resulting self- consistent 
impurity model [ll| , for a specific site, at any temperature 
T and in the limit of Jh — > oo, depends parametrically on 
fl z (the z component of the unit vector representing the 
core-spin) at that site In order to calculate the Mn 
2p core-level PES, we add to this a single 'core-hole' level 
(labeled 'c') of positive energy e c w 647.6eV at the im- 
purity site 1 1 51 ] . The core hole has an attractive Coulomb 
interaction U p d (= -6.5 eV) with both the t £' and the l V 
electrons. The resulting Hamiltonian reads : 



(d - n)a\ai + Vi(n z )(4b + tfai) 



-(Ej T + n)m + J F {m)Vt z + (U^ne - /i)tfb 



+ (e c - fi)n c + Updn c (tfb + n e ). 
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Here, the a] s create bath electrons, representing the '&' 
electrons at the other sites of the lattice , with which 
the '6' electrons of the chosen site hybridize; and which, 
for the purposes of this paper, we have approximated as 
having a discrete grid of energies {ti\. The corresponding 
hybridization parameters Vi(fl z ) are obtained from the 
We; tt z ) determined self-consistently in the DMFT [Io|, 
Jf = 2zJfS 2 (z — coordination number) and (to) 
is the magnetization. 

Both m and n c are conserved in Hcl (Eq. ((2J). Hence, 
its eigenstates can be separately calculated in the 4 sec- 
tors: (lie, n c) = (0, 0), (0, 1), (1, 0) and (1, 1), labeled here 
as (10), (FQ) (II) and (-FT), respectively. One can there- 
fore calculate a separate spectral function A cc (uj; ni,Q z ) 
in each n£ sector, and for each £l z , as : 
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Here mo, mi refer to the many-body eigenstates of Hcl 
in the n c = and n c = 1 sectors respectively for the 



(ne, Q z ) specified. The full core-hole spectrum is given 
by the weighted average 

A^^jj^^M^n^) 



where Wo(f2 z ) and Wi(fl z ) are statistical weights 
obtained as W ne (Q z )=Z(ni, n c =0, Sl z )/Z using 
Z(ni,n c ,Q z ), the constrained partition functions, 
calculated for n c =0 and for particular values of ni and 
il z ; and Z is the total partition function. 

The Boltzmann factors in the expression for A cc (oj) 
give rise to the aforementioned FED effects in PES, which 
are very weak. However, special to our model, and hence 
to manganites, are two other, unconventional, sources of 
T and x dependence which are much larger: first, the 
statistical weights Wo and W\ are dependent on x and 
T; second, as we show below, the spectra for each (ne,fl z ) 
themselves change with x and T, with redistribution of 
spectral weights over scales of eV , because the self consis- 
tent hybridization parameters Vi(Vt z ) are strongly x and 
T dependent [lfj. Hence, for simplicity, we neglect the 
FED effects in this paper, by restricting to above to just 
the ground state in the appropriate sector. 

We have calculated the spectra using two different 
methods: (1) Hcl is single-particle like in each of the 
sectors (10), (II), (F0) and (Fl), and can be exactly diag- 
onalized for any Q z . The initial states, i.e. the many-body 
ground state (GS) for each ne, are obtained by filling up 
the single-particle levels in the 10 and II sectors up to the 
chemical potential [l6j |. The final (n c = 1) many-body 
states are obtained by creating particle-hole (p-h) exci- 
tations with respect to the corresponding ground states, 
in the F0 and Fl sectors We find that the PES 

spectrum is dominated by the single p-h channel (spec- 
tral weight > 95%) which is calculable [l6| to very high 
accuracy even for a dense grid for the bath electrons. (2) 
More involved calculations including contributions from 
all p-h channels, but limited to using only 21 bath states, 
have been carried out using the Lanczos recursive algo- 
rithm [l^ . In both cases, as is standard practice, the 
discrete spectra obtained have been broadened using a 
Gaussian broadening with a ~ 0.1 eV. The single and 
all p-h spectra are practically identical, and the small 
missing weight (3-5%) in the single p-h channel is visible 
only on close inspection. All the calculations reported in 
this paper are for model parameters appropriate for the 
LBMO thin- film samples of Tanaka et al. [6J] . 

Our results of the doping variation of the PES, ob- 
tained at T = in the fully spin-polarised FM phase 
(Cl z = 1), are shown in Fig. 1. Fig.s 1(a) and 1(b) show 
the spectra separately in the two ni sectors and for both 
single and all p-h channels, for x = 0.1, 0.2, and 0.3. 
The fully averaged spectrum (Fig. 1(c)) thus has a main 
peak with two shoulders, one on each side , separated 
by about 6.5 eV (~ |t/ p d|) from a high energy correla- 
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FIG. 1: (Color online) (a),(b) Calculated core-level PES spec- 
tra at T — 0, shown in two parts for clarity, in the single p-h 
(SPH) and all p-h (APH) channels, for the n; = and m = 1 
sectors and for dopings 0.1, 0.2 and 0.3. (c) The fully aver- 
aged APH spectra for the three dopings (normalized to the 
main peak), showing only the main peak and the low energy 
shoulders. 
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FIG. 2: (Color online) (a) Magnified view of the T dependence 
of the calculated spectra averaged over the two ni sectors as 
well as the core spin angles, in the single p-h (SPH) channel 
for x=0.15. (b) The T dependence of the integrated shoulder 
intensity in the difference spectrum with respect to the para- 
magnetic spectrum, compared with M 2 (T), the square of the 
(DMFT derived) magnetization. M(T) is also shown. 



tion satellite (not shown in Fig. 1(c)) [13] ■ As seen in 
Fig. 1(c), with increasing x till about 0.3, the main peak 
intensity decreases, while those of the shoulde rs ( and of 
the correlation satellite) , increase substantially [l7fl ■ This 
is similar to what is seen in the core-level PES data in 
LSMO @, In the data on LBMO @ (where x is 

not varied), only the low energy shoulder is distinguish- 
able. We believe that this is because LBMO has a lower 
x c for the Ferro-Insulator (FI) to FM transition (0.05 
as opposed to 0.16 in LSMO), and hence [l(| a smaller 



(Ejt/Dq). In that case the small T dependent changes 
in the higher energy shoulder, separated only by ~ Ejt 
from the main peak which is rather broad in the experi- 
ment, are harder to distinguish (see Fig. 2(a) where this 
effect is simulated using a broadening of 0.3 eV). 

Fig. 2(a) shows a magnified view of T dependence of 
the fully averaged spectrum in the single p-h channel for 
cc=0.15. Clearly, spectral changes as T increases are sim- 
ilar to the changes as x decreases, precisely as observed in 
LSMO (Ij. In Fig. 2(b) we show the T dependence of the 
integrated shoulder intensity in the difference spectrum 
with respect to the paramagnetic spectrum, and find that 
it tracks M 2 (T), the square of the magnetization (as cal- 
culated using the DMFT). Such a correspondence has 
indeed been observed experimentally in LBMO [6j, and 
is not reproducible by conventional rigid-band or cluster 
calculations. 

All the major features of our calculated spectra can 
be understood from the single p-h channel contributions. 
Fig. 3 depicts the initial state and two important types of 
final states in each of the ng sectors in the metallic regime, 
with the bare bandwidth D = 1.3 eV, E JT = 0.29 eV. 
The band occupancy (n&) is small and the chemical 
potential, pinned at /i = —Ejt, lies close to the effective 



bottom edge, —D, of the l V band [la ]. 

When n( = (Fig. 3(a)) the local level in the ini- 
tial state (n c — 0) is at zero and hybridizes sparingly 
with the levels near the band edge, as the corresponding 
hybridization amplitudes are small. Thus the local l V 
character of the filled levels in the initial state is small. 
When n c — 1, the local level is pulled down by an 
amount \U p d\ — 6.5 eV, and becomes substantially oc- 
cupied in the final GS. The local l V character of the oc- 
cupied levels in the band is again very small. Hence in 
the ni = sector, the transitions involving p-h excita- 
tions with energies close to and above the spectral edge 
corresponding to the GS-to-GS transition — > [i + in 
Fig. 3(a)), at an energy ~ (U pd + e c - 2fi) S 641.68 
eV, have low intensity. The dominant contribution comes 
from the transitions to final states where the local 'b' elec- 
tron is excited to levels just above the chemical potential 
(b — > /i + in Fig. 3(a)). The corresponding edge is at an 
energy - (-E JT + e c - 2/x) ^ 647.89 eV and leads 
to the correlation satellite in Fig. 1(a). 

When m = \ (Fig. 3(b)), the local '6' level is pushed 
up, by an amount Udd = 5 eV , in the initial state (n c = 
0). This further reduces the amount of 'V mixing and 
hence the local '6' character of the occupied levels in the 
initial state. When n c = 1 (Fig. 3(b)), the '6' level is 
at (Udd + Upd) — — 1.5 eV, much closer to the band edge 
than in the ng = sector. Nevertheless, the GS when 
n c = 1 still has a substantial occupancy of the local '6' 
level. Hence, the ^~ — > ^ + transitions (Fig. 3(b)) again 
have a small intensity in the spectrum. The edge is now 
at an energy ~ (U dd + 2U pd + e c - 2p) = 640.18 
eV, 1.5 eV below the edge in the ne — sector. Just 
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FIG. 3: (Color online) Schematic depiction of the GS configu- 
rations (with the occupied levels shown shaded) in the initial 
state (n c = 0) and the ground and excited state configura- 
tions in the final state (n c = 1) for the two sectors rii = 
(a) and m = 1 (b). Transitions to final states with the local 
'6' level occupied have little intensity. The dominant transi- 
tions are to final states where the local '6' electron has been 
transferred to just above the Fermi level. 



as in the ni = sector, the main contribution to the 
spectrum comes from the b — > fi + transitions (Fig. 3(b)), 
beyond the edge at an energy [h| ~ (U p d — Ejt + e c — 
2//) = 641.39 eV, hence below the GS-to-GS edge from 
the n,£ = sector. Associated with the edge spectra are 
edge singularities [2(| and tails due to p-h excitations, 
which, when smoothed out, give rise to the shoulders 
with asymmetric 2l| lineshapes. 

As one increases x or decreases T, the '6' bandwidth 
increases 



18[. However, /it still remains close to —Ejt, so 
that the filling and the local '6' character of the occupied 
band levels near fj, increase in the initial state. Hence, one 
gets a steady transfer of spectral weight from the features 
where the local '6' level is unoccupied in the final state, 
to those where it is occupied in the final state, as seen in 
the bare or un-averaged spectra (Fig. 1(a)). The spectra 
shown in Fig. 1(b) are averaged over the contributions 
from the two ni sectors (and additionally over Vl z for 
Figs. 2(a) and 2(b)), with statistical weights which are 
themselves functions of x and T. The net effect is that 
the two shoulders on the two sides of the main peak, 
arising from smoothed out edge spectra as shown above, 
increase in intensity with increasing x or decreasing T, 
in agreement with experiments. 

In conclusion, we have presented Mn 2p core-level PES 
calculations by extending a new model for manganites 
[Io| that takes into account the simultaneous presence 
of strong electron-lattice, spin-spin and charge-charge in- 
teractions. Our results reproduce, for the first time, the 
unusual redistribution of spectral weight over several eV 
upon varying x and T, and a correspondence between the 
effect of increasing x and decreasing T, as experimentally 
observed [1, @. 
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